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Charts and formulas are presented from which bivariate normal proba- 
bilities may be computed. Formulas involving the bivariate normal are given 
for the solution of a problem of handling errors in testing and measuring. 
These formulas include, in addition to previously published cases, two new 
cases. In one, the product is not necessarily centered relative to two-sided 
specification limits; in the other, one-sided specification limits are con- 
sidered. 

I. INTRODUCTION 

Many manufactured products are 100 per cent tested with the idea of 
insuring that each unit of product meets the performance specifications. 
The general procedure is to set the test specification limits, which are the 
limits used for product unit acceptance, at or arbitrarily near the per- 
formance specification limits. These are established by engineering re- 
quirements and require a test-set accuracy and precision of a specified 
amount with respect to the performance specification limits or some 
nominal value. Eagle 1 and many others have studied the problem of 
locating test specification limits with respect to performance specifica- 
tion limits under various conditions of test-set precision when testing 
itself is subject to random error. He pointed out that, when random errors 
of testing exist, two types of errors or mistakes can occur which should 
be taken into account in setting test specification limits. The first error, 
called consumer's loss (CL), is defined as the probability that noncon- 
forming product units will be accepted. The second error, called pro- 
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ducer's loss (PL),* is denned as the probability that conforming product 
units will be rejected. 

The problem of locating test specification limits with respect to per- 
formance specification limits may be considered heuristically as in Fig. 
1 , where a and d are performance specification limits, and b and c are 
test specification limits. Consider a product unit which just fails to meet 
the lower performance specification limit and therefore is nonconforming 
and should not be accepted. The chance that a product unit of this value 
will be accepted by the test set is shown by the shaded area under the 
test-set-error distribution curve to the right of b, the lower test specifi- 
cation limit. This then is a part of the consumer's loss (CL), and the sum- 
mation of similar considerations for all product units outside the per- 
formance specification limits constitutes the CL. Consider now a product 
unit with value at c, well within the performance specification limits, 
which should be accepted. The shaded area under the test-set-error dis- 
tribution curve to the right of c, the upper test specification limit, shows 
that there is a 50 per cent chance of the test set reject ling a product unit 
of this value, and thus this is a part of the producer's loss (PL). The sum- 
mation of similar considerations for all product units inside the per- 
formance specification limits constitutes PL. Different settings of the 
test specification limits, b and c, with respect to the performance speci- 
fication limits, a and d, and/or different spreads (standard deviations) 
of the product and test-set distributions would lead to different pro- 
ducer's and consumer's losses. 



--PERFORMANCE SPECIFICATION H 




Fig. 1 — Diagram showing effect of test-set errors on test acceptance limits. 

* Consumer's loss and producer's loss are called consumer's risk and producer's 
risk by Grubbs and Coon 2 and others. The terminology used here is an effort to 
avoid confusion with present standard statistical quality control terminology. 
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Eagle presented graphs of PL and CL which assume that the product 
distribution and test-set-error distribution are both normal (Gaussian) 
and also that the mean of the product distribution is midway between 
the performance specification limits and also midway between the test 
specification limits. Thus, these two sets of limits are symmetrical with 
respect to the mean of the product distribution. He also presented for- 
mulas based on Pearson's tables 3 for computing PL and CL. Pearson's 
tables are now out of print and are not readily available to many who 
may wish to make these calculations. However, the National Bureau of 
Standards 4 plans to reissue Pearson's tables with some extensions. 

Hayes 5 and Wiesen and Clark 6 have presented additional graphs of 
PL and CL for the same conditions as those considered by Eagle. In the 
present paper, the solution is given in formula form without the product 
distribution centering requirement. Also, one-sided test and performance 
specification limits are considered. 

Grubbs and Coon 2 also considered the problem of setting the perform- 
ance and test specification limits for a centered product distribution. 
They found the locations which minimized the sum CL plus PL, and also 
the locations which minimized total cost when the consumer's loss was 
subject to a given cost and the producer's loss subject to another cost. 
In this paper, formulas are given for the same cost assumptions, but for 
both one-sided and not necessarily centered two-sided performance and 
test specification limits. 

Tingey and Merrill 7 considered the problem of minimizing the total 
cost when the cost to the consumer of accepting nonconforming product 
units varied with the degree of nonconformance. They presented a table 
of constants for constructing test specification limits under these circum- 
stances. Formulas (20) and (21) are given below for computing the total 
cost under these assumptions. 

Eagle 1 and Owen 89 showed that the bivariate normal distribution un- 
derlies the above problems. Consequently, a method for computing bi- 
variate normal probabilities is first considered here. Although this 
method is discussed relative to the present problem, the method is per- 
fectly general and can be used wherever bivariate normal probabilities 
are required. 

II. COMPUTATION' OF BIVARIATE NORMAL PROBABILITIES FROM THE CHARTS 

Let 

«* *» ■ wf^ /; r - [- \ o^s^i** « 

M(h, k; p) is the probability that a normal random variable X with 
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mean zero and variance one is greater than h, and another normal ran- 
dom variable Y with mean zero and variance one is greater than k, where 
p is the correlation between X and Y; that is, Pr(X > h, Y > k) = 
M(h, k; p). 

It is convenient also to have a functional notation for the univariate 
normal integral. To this end, define 



m = ^h L exp vrn 



Volumes of the bivariate normal over other rectangular regions may be 
expressed in terms of the M- and G-i unctions. In terms of standardized 
variables (zero means and unit variances), 



Pr(X < h, Y < k) = M(-h, -k; p) 

= G(h) + G(k) + M(h,k; P ) - 1, 
Pr(X <h,Y>k) = M(-h, k; -p) = G(-k) - M(h, k; p) 

= G(h) - M(-h, -k;p) 
Pr(X > h, Y < k) = M(h, -k; -p) = G(-h) - M(h, k; p) 

= G(k) - M(-h, -k;p). 
Ref . 8, which is an elaboration of Ref . 9, shows that 

M(h, k; P )=M (h, 0; ^ ~ *)(** *) ) 
\ Vh 2 - 2phk + W 



(2) 
(3) 

(4) 



(5) 



where the upper choice is made if hk > or if hk = but h + k ^ 0, 
and the lower choice is made otherwise; and sgn h = +1 if h ^ 0, and 
sgn h = — 1 if h < 0. Note that sgn h as used in (5) is a multiplicative 
factor; i.e., it affects the sign of the quantity it operates on but not the 
absolute value of the quantity. 

Equation (5) means that bivariate normal probabilities with any 
limits can be computed from a table of the bivariate normal probabilities 
where one of the limits is zero. Figs. 2 through 5 are graphs of these 
probabilities, i.e., of M(h, 0; p). Figs. 4 and 5 can be obtained by rotat- 
ing Figs. 2 and 3, respectively, 180 degrees and replacing h by — h and 
p by —p. However, to avoid confusion both figures are included. The fol- 
lowing examples illustrate the use of these graphs. 
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Example 1: The following probabilities may be read from Figs. 2 
through 5 : 

M(l, 0;0.5) = 0.127, 
il/(l, 0; -0.5) = 0.031, 
iV(-l,0;0.5) = 0.469, 

If (-1,0; -0.5) = 0.373. 
Example 2: Find il/(1.96, 1; 0.56). 

Solution: 

\A 2 - 2phk + k* = V2T6464 = 1.6268. 

il/(1.96, 1;0.56) = M(l. 96, 0; 0.060) + M (1, 0; -0.861) 

= 0.014 + 0.002 

= 0.016. 

The value computed by more elaborate methods 8,9 is 0.016. 
Example 8: Find M(0A, 0.1; -0.5). 
Solution: 

\/h 2 - 2phk + k 2 = V21- 

il/(0.4, 0.1; -0.5) = il/(0.4, 0; -0.655) + M (0.1, 0; -0.982) 

= 0.069 + 0.015 

= 0.084. 

The value read from Pearson's tables 3 is 0.0836. 
Example .',: Kind .1/(1, -2; 0.7). 
Solution: 

Vh? - 2phk +T 2 = 2.793. 

.1/(1, -2; 0.7) = iU(l, 0; 0.967) + M(-2, 0; 0.859) - \ 
= 0.159 + 0.500 - 0.500 
= 0.159. 

The value computed from Pearson's tables 3 is 0.1587. 
Example 5: Find Af(l, -2; -0.5). 
Solution: 

y/K- - 2 P hk + A 2 = 1.732. 

M(l, -2; -0.5) = ilf(l, 0; 0.866) + M(-2, 0; 0) - h 

= 0.157 + 0.489 - 0.500 

= 0.146. 

The value computed from Pearson's tables 3 is 0.1454. 
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III. ERRORS IN TESTING AND MEASURING 

Before giving the formulas for the problem of errors in testing and 
measuring, let the following quantities be defined. 

Let P = a random variable (normally distributed with mean n and 
standard deviation <r p ) which describes the true product values. 

Let T = a random error of measurement (test-set error) normally 
distributed with mean X and standard deviation <r t . 

The term X is called the bias of the measuring instrument. In Refs. 
1, 2, 5 and 6, X was taken to be equal to zero. This is not necessary, as 
only a slight change in the definition of the test specification limits 
takes care of this. (See the definition of 61 and b 2 below.) Of course, in 
most cases where X is known, the test instrument would be recalibrated 
and the bias eliminated. There may be situations, however, where this 
would not be practical and it really is no problem to carry X as an extra 
parameter. Also, in assessing the problem of an unknown bias, it is con- 
venient to have the parameter X available (see Example 8 below). 

Let S = the observed measurement = P + T, so that, if T and P 
are taken to be independent, S is normal with mean n + X and variance 
c s 2 = <Jp + o~t 2 - Also, P and S are correlated, and it will be shown that 
the correlation is equal to a p /a s . 

The expectation or mean of a random variable X is indicated sym- 
bolically by E(X). The correlation between two variables X and F is 
defined as 

<t x <Ty 
For P and S, then, the numerator of the correlation is 

E[(P - m) (S - n - X)] 

= E[(P - M ) (P + T - M - X)] 

= E[P 2 + PT - Pn - P\ - ftP - fiT + /x 2 + /xX]. 

Now, E(X 2 ) = <r A - 2 + m 2 ; E (constant) = constant; E(ZF) = E(X)E(F), 
if X and Fare independent; E(X + Y) = E(X) + E(F) and E(X) = fx x . 
The correlation between P and S is then 

o-p" + n~ + ^X — n — fx\ — (x — juX + m + mX _ o P 
Pps — • 

<7 p O\, <T S 

Let &i and kz be defined so that the performance specification limits 
are at n -f- k\<r p and n — kyj p . For a lower performance specification 
limit only, take h\ = 00 ; and for an upper performance specification 
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limit only, take fc 2 = °° . In Refs. 1, 2, 5 and 6, the fe's were taken to be 
equal which meant the product distribution was assumed to be centered 
with respect to the performance specification limits. 

Let b v and b 2 be defined so that, when the test specification limits are 
placed at /i + X + h<r P — ^t and at n + X — k 2 <r p + b 2 a t , either a 
desired consumer's loss or a desired producer's loss is not exceeded, or 
some combination of these losses is not exceeded. 

In this paper, n,\,<r p and <r, are assumed to be known. When any of 
them are unknown, experiments must be run to establish their values. 
In this connection, the reader may wish to read a paper by Grubbs, 10 
which gives methods for estimating <r p and c t when the variance of each 
observation is a linear function of o v ~ and <sc. Fig. 6 shows diagram- 
matically the relationship between p, X, tr p , <r t , &i , b% , ki and k 2 . 

fj, = MEAN OF PRODUCT DISTRIBUTION 
<Tp = STANDARD DEVIATION OF PRODUCT DISTRIBUTION 
0\ = STANDARD DEVIATION OF TEST- SET ERROR DISTRIBUTION 

\ = TEST- SET BIAS 



|^. PERFORMANCE SPECIFICATION 

I 

k---k 2 0p— ■»+*■ k, <Tp 

M TEST SPECIFICATION »j 

I ,-\+k 2 <r p -D 2 cr t | 

I I IT 



I 

"1 




/i-k 2 o- p — - 
(j.+\- k 2 £Tp + b 2 cr t ' 



//. + k,(r p -b,<r t -'' 
H.+ A.+ k,Op-b 1 o- t -' 



'-/i+k,(rp 



Fig. 6 — Diagram showing various constants defined for errors in the testing- 
and-measuring problem. 
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From the definition of consumer's loss (CL) it is clear that 

CL = Pr[P > n + h<r p or P < /x - k 2 a p 

(6) 
and n + X - (k 2 a p - b 2 <r t ) < S < y. + X + {k x a p - b&t)]. 

Since the correlation between P and S is given by 

P P s = p = <t p /(Ts = o-p/s/o-p 2 + (Jt 2 , 

and P and S are assumed to be bivariate normal, the above probability 
can be expressed in terms of the M-i unction. The first step, as with the 
univariate normal, is to standardize the variables, i.e., subtract the mean 
and divide by the standard deviation. In order to save space and give 
compact formulas, additional notation is introduced. Let 



ki<T p — blCTt , kiffp — b-2<T t 

* = 77 2 _t_ 2 aild «* = —J 2 , , ■ 

V <Jp + <Tt 2 VOp + <T t 2 



Then, 



C 



(L = Pr ( - £ > fa or £ < - k 2 



j ^ & ~ M — X ^ 

and — g 2 < < q\ 



- Pr (£^-£ > fa and g - M " X > - q 2 ) 

\ <T P (T s / 

— Pr I > fa and > (h ) 

\ (Tp <7 a / 

+ Pr (£j=J? < - fa and g ~ M " X > - q. 



- Pr (*LZ± < - fa and g - M = X > q\ 

Now (2), (3) and (5) are used to reduce CL to its most easily used form : 

CL = M(fa, -q 2 ; P ) - M(h , q, ; p) 

(7) 
-M{k 2 ,q 2 ;p) +M(h i} - ffl ;p). 

The producer's loss (PL) is given by 

PL = Pr Lw — &2<r p < P < ju + fcio-p and 

(8) 
S < p. + X — (fao- p — &2<r<) or £ > /x + X + (ki<r p — &i<r«)]. 
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If the same procedure as that with consumer's loss is followed, it can be 
shown that (8) reduces to 

PL = CL + G(h) + G(k 2 ) - G(qi) - G(q 2 ). (9) 

For one-sided specification limits the above formulas reduce to: 
For an upper limit only: 

CL = Pr(P > n + ki<T„ and S < y. + X + h<r p - bia t )', 

hence 

CL = Gi-h) - M(k 1 ,q 1 ; P ) (10) 

and 

PL = Pr(P < m + hffp and S > n + X + k x a p - 6un); 

and hence, 

PL = G(- qi ) -M(k ltqi ; P ). (11) 

For a lower limit only : 

CL = Pr(P < n - k 2 <r P and S > n + X - /c 2 o- P + b 2 <T t ); 

hence 

CL = G(-kz) - M(k 2 , « 2 ; p) (12) 

and 

PL = Pr(P > n — fc 2 <T P and »S < n + X — kxr p -f 620- <); 

hence 

PL = G(-g 2 ) - M(k 2 ,q 2 ; P ). (13) 

The following examples illustrate the use of (7) through (13). 

Example 6: Suppose the performance specification limits are at n + 2<r p 
and at n — 3<r p . Suppose also that the test specification limits have been 
set at n + X + 2<r p and at m + X — 3tr p + a t . If <r,/<r p = 0.5, deter- 
mine the producer's and consumer's losses of this procedure. Since the 
bias is known in this example and has been allowed for in setting the 
test specification limits, X does not appear in the calculations of PL and 
CL. (See example 8 for the case where X is not known.) 

Solution: Here, k x = 2, k 2 = 3, 61 = 0, and b 2 = 1. Then, 

ha p - hat 2-0 

91 = V77T7? = vns = L789 ' 
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k 2 <r p - bi<r t 3 - 0.5 
92 = Vctf + rf = VL25 " 2 - 236 



and 



p = ; Cp . = -7= = 0.894. 



W + <r, 2 " VL25 
From (7), 
CL = M(2, -2.236; 0.894) - i¥(2, 1.789; 0.894) 

-M(3, 2.236; 0.894) + M(S, -1.789; 0.894). 
Application of (5) results in 
CL = M(2, 0; 0.976) + M (-2.236, 0; 0.970) - 0.5 - il/(2, 0; 0) 

- M (1.789,0; -0.447) - M(3,0; 0.316) - M (2.236,0; - 0.707) 
+ M(3, 0; 0.958) + M(- 1.789, 0; 0.985) - 0.5 

= 0.023 + 0.500 - 0.500 - 0.017 - 0.004 - 0.000 

- 0.000 + 0.000 + 0.500 - 0.500 
= 0.008. 

From (7), 

PL = 0.008 + G(2) + G(3) - G( 1.789) - G(2.236) 
= 0.008 + 0.977 + 0.999 - 0.963 - 0.987 = 0.034. 

Example 7: The circumference of a product has a mean value of 28.5 
inches with a standard deviation of 0.5 inch. Only circumferences less 
than 29 inches are acceptable. The device for measuring the circum- 
ference is known to be biased so that on the average it measures 0.1 
inch too small with a standard deviation of 0.2 inch. If the upper test 
specification limit is set at 29 inches, and there is no lower limit, what 
are the producer's and consumer's losses? 

Solution: Here, n = 28.5, a p = 0.5, X = —0.1 and at = 0.2. 
Hence, for the performance specification limits, 

M + fciffp = 29, 

or 

28.5 + fci(0.5) = 29. 
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Solving this equation results in A'i = 1. For the test specification limits, 

H + X + h<r p - bun = 29, 
or 

28.5 - 0.1 + 1(0.5) - foi(0.2) = 29. 
Solving this equation results in 6 X = —0.5. Hence 



ki<r p - bur i 0.5 + 0.01 
Qi - 

and 



9l ~V^+tf- Vow ~ 1114 ' 



0.5 

= 0.928. 



P V<r P 2 + it- V0.29 

From (10), 

CL = G(-l) - 3/(1, 1.114; 0.928) 

- 0.159 - il/(l, 0; - 0.447) - #(1.114, 0; 0.08.3) 

= 0.159 - 0.036 - 0.073 = 0.050. 

From (1 1 ), 

PL = G(- 1.114) - 1/(1, 1.114; 0.928) = 0.133 - 0.109 = 0.024. 

Thus, for the measuring device considered here, the consumer's loss is 
0.050 and the producer's loss is 0.024. 

Example. 8: If the bias X is assumed to be zero when it is actually posi- 
tive, the effect is to accept more products in the lower tail of the product 
distribution and to reject more in the upper tail. The effect on CL and 
PL is the same as if kia p had been decreased by X and A' 2 o- p increased by 
X, i.e., as if both test specification limits had been decreased by the quan- 
tity X. 

As an example, the performance specification limits for a product are 
at m + A-jo-p = 90 and /z — A- 2 o> = 80. Also, p = 85, at = 1 and a p = 2. 
Then fej = k« = 2.5, and Table I of Grubbs and Coon 2 can be used to 
find bi = b 2 so that the consumer and the producer accept equal losses. 
The table gives &i = 62 = —0.5902 and results in a consumer's loss and 
producer's loss each equal to 0.0061 or a total loss of 0.0122. Now suppose 
that the tester is biased one unit high, i.e., X = 1. What effect does this 
have on the producer's loss and the consumer's loss? 



568 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1959 

Solution: The test limits were set at 

M + A-io- P - &*r, = 90.5902 
and 

H - kocTp + b 2 <j t = 79.4098, 

whereas, if the bias had been taken into account, they should have been 
set at 

M + X + k l( r p - hurt = 91.5902 

and 

M + X - k 2 cT p + b- 2 <T, = 80.4098 

to maintain the equality of PL and CL. To find the effect of this error, 
it is only necessary to decrease ki<r p by X; i.e., redefine A'i so that 
ki<r p = 5 — X = 4 and hence take fa = 2, and redefine /c 2 so that 
k7,<T p = 5 + X = 6 and hence take k 2 = 3. Now the formulas for CL 
and PL [(7) and (9)] are used with the following parameters: 

k x = 2, fc 2 = 3, 

hi = -0.5902, b 2 - -0.5902, 

ff p =2, <T| = 1. 

CL = il/(2, -2.947; 0.894) - M(2, 2.053; 0.894) - M(3, 2.947; 0.894) 

+ M(3, -2.053; 0.894) 
= il/(2, 0; 0.982) + il/( -2.947,0; 0.962) - 0.500 - M(2,0; -0.281) 

- M (2.053,0; -0.180) - il/(3,0; -0.193)- il/(2.947,0; -0.267) 

+ Af(3, 0; 0.962) + M (-2.053, 0; 0.982) - 0.500 
= 0.023 + 0.500 - 0.500 - 0.006 - 0.007 - 0.000 - 0.000 + 0.001 

+ 0.500 - 0.500 

= 0.011. 

PL = 0.011 + (7(2) + <?(3) - (7(2.053) - (7(2.947) 

- 0.01 1 + 0.977 + 0.999 - 0.980 - 0.998 = 0.009. 

Hence, because of the bias in the test set, the consumer's loss is increased 
from 0.006 to 0.011, and the producer's loss is increased from 0.006 to 
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0.009. If the same test specification limits are used over a period of time, 
it would probably be worthwhile to compute the consumer's loss and 
producer's loss for a whole range of biases. These could be plotted on a 
graph showing the assumed bias on the abscissa and the producer's or 
consumer's loss on the ordinate. 

IV. SPECIAL CRITERIA FOR DETERMINING SPECIFICATION LIMITS 

Following Grubbs and Coon, 2 if the condition is that the producer and 
consumer accept the same or equal losses in rejecting a conforming prod- 
uct unit and in accepting a nonconforming product unit, the solution is 
obtained by setting CL = PL in (9) for two-sided test and performance 
specification limits, or setting CL = PL in (10) and (11) or (12) and (13) 
for one-sided test and performance specification limits. This results in 
</i = A'i and q 2 = k 2 or 

, _ fcio-p — fci \/<t p 2 + a? / 14 s 



at 



and 



b 2 = fog>-*» vV + g''. (15) 



at 



Equations (14) and (15) may be solved for b x and b 2 , respectively; or, 
if r is set equal to <r p /<rt , the value of b satisfying (14) and (15) may be 
read from Table I (p. 16) of Ref. 2. The values for consumer's loss and 
producer's loss, however, will have to be calculated from (7), (9), (10), 
(11), (12) and (13). 

Another criterion discussed by Grubbs and Coon is to assume that the 
sum of the consumer's and producer's losses is to be a minimum. Then 
the b's for two-sided test specification limits are obtained by solving the 
equations 



and 



/ - k 2 - P qA /- h + Mi \ 

r / - fei - pg 2 \ r ( - h + pg2 \ 



Again following Grubbs and Coon, if fci ^ 1.5, k 2 ^ 1.5, <r p g£ <r< and 
the b's are negative or small positive, then the first integral in each of 
the above equations is nearly zero, so that approximate solutions to the 
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above equations are fcj = pqi and £ 2 = pqi, or 

i,a-^ (16) 

and 

feS*-^- 1 . (17) 

The /> values can be read from Table II (p. 16) of Ref. 2, but the con- 
sumer's and producer's losses have to be computed from the formulas 
given in this paper. Equations (16) and (17) for the b's are exact for one- 
sided test specification limits. 

Let Cci be the cost of accepting a nonconforming product unit, and 
C p i be the cost of rejecting a conforming product unit. Then the values 
for 6i and 6 2 that minimize the total cost are those that satisfy the equa- 
tions 



and 



f -h- pq{ \ r /-fe i + PQ i\ _ C p i 

9 \ VT=?) + WW / 5 + c, 

/ -/>-i ~ pq* \ (~ h + pq A _ C pt 



G 
and 

G 



Again, if fc x ^ 1.5, fc 2 ^ 1.5, <r p ^ <r t and the 6's are negative or small 
positive, then the first integral in each of the above equations is nearly 
zero and the b's may be obtained approximately from the equations 

/ -hat — bia p \ ^ C p i , lg N 

\ V Op 2 + <*? ) Cd + C p i 

( "/"'""? ) ^ ^r^V ■ (19) 

\ VOp + OT / Ce| -+- tpi 

In this case the b's may be obtained approximately (for two-sided test 
specification limits) by setting the limit of the integral to that devi- 
ate of a univariate normal which corresponds to the fraction 

Cpi 
Cd -\~ C p i 

Tingey and Merrill 7 also consider this problem, and they allow the cost 
of accepting a nonconforming item to be different for the two tails of the 
product distribution. They give a short table of values of 6i and b 2 (their 
Table II) under these conditions. Note that the subscripts 1 and 2 must 
be interchanged to enter their table with the formulas given here. 
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For a one-sided test specification limit the b is obtained exactly by 
(18) and (19). The corresponding producer's and consumer's losses may 
then be calculated from (10), (11), (12) and (13) given in this paper. 

As mentioned in the introduction, Tingey and Merrill 7 consider the 
case where the producer's loss is constant but the consumer's loss varies 
with the degree of nonconformance. They give a table of values of 6i 
and b 2 which will minimize the total cost under these conditions. Note 
that they use the subscript 1 for the lower tail of the product distribu- 
tion and 2 for the upper tail, while the reverse has been used here. Hence, 
if the formulas in this paper are being used, interchange 6i with 6 2 and 
A'i with A- 2 when entering Tingey and Merrill's table. They do not give 
values for the total loss, but their formula for the total loss may be re- 
duced to a computing form as follows: 

Let Cd" = cost to the consumer of accepting a nonconforming prod- 
uct unit from the upper tail only of the product distribution. Then, if 
the cost to the consumer varies with the degree of nonconformance, this 
is defined mathematically to be 

(~i u _ 1 P- A'KTp ~u 

where C" is the unit cost associated with the acceptance of nonconform- 
ing product in the upper tail of the product distribution. Then, in the 
notation used here, Tingey and Merrill's formula can be reduced to 

^ = -M/(/u , -q-i ; p) + hM{h , Qi ; p) 



P 



exp(-p-f)G 



— Q\ + p/>'i 



+ 






Hence, with the aid of Figs. 2, 3, 4 and 5 and a table of the univariate 
normal distribution function and its derivative, it is possible to compute 
C c ". A similar formula for the consumer's cost for the lower tail of the 
product distribution may be obtained by interchanging Ai with fcs and 
interchanging </i with q* . 

If the cost per product unit to the producer of rejecting conforming 



+ 
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product is constant, the total cost to the producer may be found by mul- 
tiplying the PL obtained from (9) [where the CL in (9) is computed from 
(7)] by the per product unit cost. Then the total cost to both producer 
and consumer is obtained by adding (a) the consumer's cost for the up- 
per tail of the product distribution, (b) the consumer's cost for the 
lower tail of the product distribution and (c) the producer's cost. It is 
this total cost which is minimized by the values of &i and b 2 found in 
Tingey and Merrill's Table I. 

Consider now the case where the producer's cost of rejecting a con- 
forming product unit is proportional to the degree of nonconformance, 
but only the upper tail of the product distribution is to be considered. 
The consumer's cost can be obtained by putting q 2 equal to infinity in 
(20). The result is 

<gL = -kiG(-h) + hM(h ,qi; P ) 

-<m exp{ -* qi)G \vr^) (21) 

As before, a similar formula holds if a lower tail only is considered and 
may be obtained by replacing fci with k 2 and q\ with q 2 . 
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